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I. INTRODUCTION 



Gravitational Radiation from Collapsing Magnetized Dust. II 
— Polar Parity Perturbation — 

Hajime Sotanj^ 

Theoretical Astrophysics, University of Tubingen, 
Auf der Morgenstelle 10, Tubingen 72076, Germany 
(Dated: April 9, 2009) 

With the gauge-invariant perturbation theory, we study the effects of the stellar magnetic fields 
on the polar gravitational waves emitted during the homogeneous dust collapse. We found that the 
emitted energy in gravitational waves depends strongly not only on the initial stellar radius but 
also on the rasio between the poloidal and toroidal magnetic components. The polar gravitational 
wave output of such a collapse can be easily up to a few order of magnitude larger than what we 
get from the nonmagnetized collapse. The changes due to the existence of magnetic field could be 
helpful to extract some information of inner magnetic profiles of progenitor from the detection of the 
. gravitational waves radiated during the black hole formation, which results from the stellar collapse. 

^ ' PACS numbers: 04.25.Nx, 04.30.Db, 04.40.Dg 
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' The direct observation of gravitational waves is a significant attempt in theoretical and experimental physics. In 

^h' order to accomplish this great purpose, several ground-based laser interferometric detectors with kilometer-size arm, 
such as LIGO, TAMA300, GEO600, and VIRGO, are currently operating and the next-generation detectors are also 
on the menu In addition to the ground-based detectors, the projects to launch a detector in space, such as LISA 
01 and DECIGO 0], are in progress. The importance to detect gravitational waves directly is that one can obtain 
a new method to see the Universe by gravitational waves, which is called "gravitational-wave astronomy". In fact, 
by using the gravitational waves associated with the oscillations of coinpact object, it is possible to determine the 
stellar radius, mass, equation of state (EOS), and so on 

[1 S, i, 0, B @, 0, El- Additionally, with collecting the 
\^ ' observational data of gravitational waves, we may able to verify the gravitational theory, and new physics at a high 

; density or high energy region. 
T-H For both ground-based and space detectors, the nonspherical stellar collapse is one of the most promising sources 

.,— ^ ; of gravitational waves. With their sensitivity, the black hole formation with stellar mass is a target for the ground- 
' based detectors, while the space detectors might detect signals radiated from the creation of supermassive black holes 
0^ ! [i-2, 13, 14]. The most hopeful approach to calculate these gravitational waves emitted during the black hole formation, 
' is the numerical relativity, i.e., via direct numerical integration of the exact Einstein and hydrodynamic equations, 
^ ■ which are fuU-nonlinearly coupled with each other. In the last decade, the numerical relativity made dramatic 
developments and it becomes possible to treat many complicated matter and spacetime [l^ [l^ [l3, [3- Still, the 
/\ computation of gravitational waves with high accuracy is not technically easy, because the gravitational waves emitted 
from the stellar collapse are very weak and sometimes they could contain unphysical noises due to gauge modes and/or 
numerical error. So in this paper, as an alternative approach, we consider the linear perturbation theory. This is 
possible to extract the weak gravitational waves with precision and also becomes a cross-check for the numerical 
results with numerical relativity. 

With respect to the calculation of gravitational waves radiated from the stellar collapse to the black hole, initially, 
Cunningham, Price and Moncrief derived the perturbation equations on the Oppenheimer-Snyder solution, which 
describes a homogeneous dust collapse jlOi], and calculated radiated gravitational waves [l^. Subsequently, Seidel 
and co-workers studied the gravitational waves emitted from the stellar collapse when a neutron star would be born 
[2l| , where they used the gauge- invariant perturbation formalism on the spherically symmetric spacetime formulated 
by Gerlach and Sengupta 22]. Iguchi, Nakao and Harada investigated nonspherical perturbations of a collapsing 
inhomogeneous dust ball [23|, which is described by the Lemaitre-Tolman-Bondi solution Further, Harada, 

Iguchi and Shibata calculated the axial gravitational waves emitted from the collapse of a supermassive star to a 
black hole by employing the covariant gauge-invariant formalism on the spherically symmetric spacetime and the 
coordinate-independent matching conditions at stellar surface, which is devised by Gundlach and Martfn-Garcfa 
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[201 . Recently, with the same formahsm, Sotani, Yoshida and Kokkotas considered the magnetic effect on the axial 
gravitational waves emitted from the collapse of homogeneous dust sphere [l^l (hereafter, we refer to this article as 
Paper I). 

In spite of many investigations of gravitational radiation from the stellar collapse with linear perturbation analysis as 
mentioned the above, they have not included the effect of the magnetic fields on the emitted gravitational waves, except 
for Paper I. It should notice that Cunningham, Price and Moncrief also dealt with the electromagnetic perturbations 
on the Oppenheimer-Snyder solution, but they did not consider the direct coupling between fluid and magnetic field 
[20| . While, the importance of magnetic effects on the evolution of compact objects has been recently realized due 
to the appearance of new instruments with high performance. One of the most remarkable examples is the discovery 
of magnetars, which are neutron stars with strong magnetic field such as B > 10^^ Gauss. With magnetar models, 
it is successful to explain the observed specific frequencies of quasi-periodic oscillations (QPOs) in the decaying tail 
of the giant flares [2^, [H, [s^]- Since there exists a quite strong magnetic field in some neutron stars, which would 
be produced after the stellar collapse, it is natural to take into account its effect on stellar collapse. And it could be 
also probable that the magnetic fields of the collapsing object are amplified during the collapse due to the magnetic 
flux conservation and they would affect the emitted gravitational waves, even if the initial magnetic field is weak. 
Actually, Paper I shows the possibility that the magnetic fields affect on the axial gravitational waves emitted during 
the dust collapse. Additionally, there exists another example showing the importance of magnetic effects in the 
evolution of compact objects, which is related to gamma ray bursts (GRBs), i.e., the short-duration GRBs could arise 
from hypcrgiant flares of magnetars associated with the soft gamma repeaters (SGRs) [3l| or collapse of magnetized 
hypermassive neutron star |32l | . 

Indeed, all examples mentioned above suggest that the magnetic fields play an important role in the stellar collapse 
and its effect should not be negligible. So in this paper, in order to explore the effects of magnetic fields on the 
emitted gravitational waves for the black hole formation, we consider the polar gravitational waves radiated during 
the collapse of homogeneous dust ball with weak magnetic field. In particular, along with Paper I, we focus on the only 
quadrupole gravitation waves, which could be more important in the astrophysical point of view. The weak magnetic 
fields are treated as small perturbations on the Oppenheimer-Snyder solution and we make an investigation with the 
covariant gauge-invariant formalism on the spherically symmetric spacetime and the coordinate- independent matching 
conditions at stellar surface proposed by Gundlach and Martfn-Garda [1^ . It should be emphasized that so far there 
is no calculation of polar gravitational waves on the dynamical background spacetime with the same method, i.e., this 
paper is first calculation. The reason why the polar gravitational waves could not be solved, might be the difficulty 
to deal with the boundary conditions at the stellar surface due to the existence of many perturbative variables, in 
contrast to the axial gravitational waves. Furthermore, as shown in the main text for ordering of perturbations, 
we consider the first order perturbations for the metric and fluid motion while the second order perturbations for 
magnetic fields, to see the magnetic effects on the emitted gravitational waves. 

This paper is organized as follows. In Sec. [Til we briefly describe the gauge-invariant perturbation theory on the 
spherical symmetric background, the background solution that we adopt in this paper, and how to introduce the 
magnetic fields as the perturbations. Next, in Sec. IIIIl we derive the perturbation equations for polar gravitational 
waves emitted during the collapse of magnetized dust sphere. Then the details of numerical procedure are shown in 
Sec. IIVI and we devote Sec. [Vjto description of the code tests. In Sec. IVI[ we show the numerical results related to 
the influence of existence of magnetic fields on the gravitational waves radiated during the formation of black hole. 
Finally we make a conclusion in Sec. IVIII In this paper, we adopt the unit of c = G = 1, where c and G denote the 
speed of light and the gravitational constant, respectively, and the metric signature is (— , +, +, +). 



II. BASIC PROPERTIES 



As with Paper I, we deal with the electromagnetic fields as the small perturbations on the dust sphere, since the 
magnetic energy is much smaller that the gravitational binding energy even if the source of gravitational waves would 
involve a strong magnetic field like a magnetar. Thus the background metric and four- velocity of fluid are 
determined as solutions of a collapsing spherical dust sphere without the electromagnetic fields. Now it is convenient 
to introduce two small dimcnsionless parameters related to strength of the magnetic field and to amplitude of the 
gravitational waves, i.e., 77 ^ \B / (GM'^ Rj'^y^^\ and e ~ |<5g^i/|, where Rs is stellar radius and we assume that the 
fiuid perturbations are also small, i.e., |(5u'^| ^ e and \Sp\ ^ e. Then the leading terms for the perturbations of t^/^J"^ 

and t^^^^ are St'^J^^ 0{e) and i5t^^*^^ ^ 0{rf), where i^t^^ and t^^*^^ express the energy-momentum tensors for the 
fluid and for the electromagnetic fields, respectively. In this paper since we focus on the effect of magnetic field upon 
the emitted gravitational waves during the stellar collapse, we omit terms of higher order such as 0{e'^) and 0{e^Tf). 
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Further with assumption that e ^ , the perturbed Einstein equations of order e are reduced to the following form 

= 87r{4^^) + <5t(f *^)} + 0{e') = 8nSt,^ + 0{e') . (2.1) 

Notice in this approximation the gravitational perturbations are driven both by the magnetic field and the fluid 
motions of the collapsing dust sphere. 

A. Gauge-Invariant Perturbation Theory 

For spherically symmetric background spacetime, the first order gauge-invariant perturbation theory has been for- 
mulated by Gerlach and Scngupta ^22^ and further developed by Gundlach and Marti'n-Garcfa 26]. In this subsection 
we briefly describe this formalism for the polar parity perturbations. 

1. Background Spacetime 

The background spacetime, which is a spherically symmetric four dimensional spacetime M., can be described as a 
product of the form Ai — M.' x 5^, where Al^ is a 2-dimensional (1+1) reduced spacetime and S' a 2-dimensional 
spheres. In other words, the metric g^jy and the stress-energy tensor tf^i, on M. can be written in the form 

g^j. = diag(gAB, r^7ah), (2.2) 
EE diag(tAB, Qr'lab), (2.3) 

where gAB is an arbitrary (1 + 1) Lorentzian metric on Ai', r a scalar on Q some function on and jab is the 
unit curvature metric on S'. Note that if the background spacetime is spherically symmetric then Q = f^^/'^ ■ Here 
and henceforth the Greek indices denote the spacetime components, the capital Latin indices the Ai' components, 
and the small Latin indices are used to denote the S' components. Furthermore, the covariant derivatives on Ai, 
Ai'^, and S' are represented by and :a, respectively. Finally, the totally antisymmetric covariant unit tensor on 

Al^ is denoted as Eab and on S' as Eat- 



2. Nonradial Perturbations 



As mentioned before, in this paper, we consider axisymmetric polar parity perturbations both for the metric 5g^^ 
and the matter perturbations St^u, which are given by 

g ^ ( MabY'^ ^tfv'"^ \ , . 



where Z^^" = r'™^ + 1(1 + l)y'™7afc/2 and y'™ stands for the spherical harmonics. With Hab, h>l\ K, G, AtAB, 
At^^\ At'^, and Ai^, the gauge-invariant variables for the nonradial perturbations are defined as 



kAB 
k 

Tab 


= hAB - {pa\b +Pb\a)^ 

= K- 2« V, 

A c c c 

= i\tAB — tAB\CP — iACP \B ~ ^BCP |A' 


(2.6) 
(2.7) 
(2.8) 


Ta 


2 

= At'^^^~tAcP^-^-QG\A, 


(2.9) 




= At' - r'QG, 


(2.10) 


^3 


EE At^ - (Q|c + 2Q«c)p^ + ^^^gC, 


(2.11) 
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where va = r^^/r and pA = h'^J^ — r'^G\A/'2 [111- Note that Ta is defined for Z > 1 and for I >2. In terms of the 
gauge- invariant variables, the Unearized Einstein equations for the polar parity perturbations are given by [2^; 



Gr 



- kcA\B — kcB\A 
+ 9AB {^V 
+ 2 {vAk\B + VBk\A + k\A\B) - 9AB 



2k 



kAB - '^gABv'^ {kED\C - kcE\D " kcD\E) 9^^ 

n 



k ^ 



-WttTab: 



k\A fcyiJ 



ka^, i - VAkg^ = -IGttTa, 



"B \A 



2v^k^A + G^^k 



\a\b 



"AB 



2v'^k 



\B 



AB 



2(v^\^ +V^V^]kAB 



k -4 \B:„,Au B A ^(^ + 1) , . 

'^A \B '^B \A + '^'^A A 



^ 167rr^ 



where TZ is the Gaussian curvature on Ai'^, and Gab and G^ are defined as 

Gab = -2 {va\b + vavb) + gABVo, 



G " = 2 



VAV^-n] , 



(2.12) 
(2.13) 

(2.14) 
(2.15) 

(2.16) 
(2.17) 

U = U^WA, 



where Vq is defined as = 2 + W' - ijU + vW^ + 3 (VK^ - U^) - r'^ and n 

W = n"^VA, F = u^F\A, and F' = n^F\A [s^l- Now if the symmetric tensor kAB is decomposed with a coordinate- 
independent way into three scalars, such as 



kAB = q{-UAUB + nAns) + <j){uAUB + nAns) + 4:{uAnB + uaub), 



(2.18) 



and to eliminate we introduce the new variable C defined as C = — fc + q, then from Eqs. p.l2p - p.lSp we can 
get the perturbation equations for the variables of metric perturbations as 



-C + C" + 2(m-c/)^'-5c. 

-k + clk" -2ciUi}' = Sk, 
q = -SttT^ 



(2.19) 
(2.20) 
(2.21) 
(2.22) 



where the source terms S"^, Sk, and are given in Appendix [XI which are without the variables for matter pertur- 
bations. Namely, we can calculate the evolutions for metric perturbations independent of the matter perturbations. 
On the other hand, the variables associated with the matter perturbations are given by 

-Sttu'^u^Tab = (fc)' + C^, (2.23) 
Sttu'^u^Tab = -k" + 2U%Ij' + G^, (2.24) 
~167rM^r4 = V' + Ca- (2.25) 

Notice that the right hand sides of Eqs. (I2.23P - (|2.25p are produced with only variables for metric perturbations, i.e., 
the matter perturbations can be determined after the calculation for the metric perturbations. The concrete forms of 
G'Y , Gu , and Gq- are described in Appendix lAl 



B. Oppenheimer-Snyder Solution 



We briefly describe the adopted background spacetime which will be later endowed with a magnetic field. We 
consider perturbations around a homogeneous spherically symmetric dust collapse described by the Oppenheimer- 
Snyder (OS) solution, whose line element inside the dust sphere is given by 

ds^ = -dr^ + R^{T)[dx^ 4- sin^ x{dO'^ + sin^ Bdcf)'^)] , (2.26) 
= R^{j])[-dif + dx^ + sin^ x{de^ + sin^ edcj)^)] , (2.27) 
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where x is a radial coordinate defined in the range of < x < xo < 7r/2 and xo corresponds to the stellar surface. 
Additionally R{r]) and r(r?) are the scale factor and the proper time of a comoving observer with the fluid, respectively, 
which are defined in terms of the conformal time r] as follows 



M , 
— o — (1 + cos 77) , 
sin Xo 

M . 
—3 — (r? + smr?), 

sm Xo 



(2.28) 
(2.29) 



where M is the total gravitational mass of the dust sphere. The energy-momentum tensor for the dust fluid is written 
as 



where p is the rest mass density given by 



pill) 



3 sin^ Xo 



(1 + cos 77) 



47rM2 

and denotes the four- velocity of the dust, described in terms of comoving coordinates as 

w'' = (5'; or u^ = R{r])-^5% 



(2.30) 



(2.31) 



(2.32) 



where 5''^ means the Kronecker delta. Also with the four- velocity u'', the spacelike unit vector defined as n"^ = —eabu^ 
is given by 



(2.33) 



where n"^ is a normal vector on the sphere whose radius is constant. Thus we have p, = U = dr^R/R?, v = 0, and 
W = cosx/{Rsinx)- Furthermore the frame derivatives arc F = dj^F/R and F' = d^F/R. The spacetime outside 
the dust sphere is described by the Schwarzschild metric, i.e.. 



ds' = -f{r)dt' + f{r)-'dr' + r^dO' + sin^ Odcj)'), 



(2.34) 



where /(r) = 1 - 2M/r. Note that with the Schwarzschild metric = (l/v7,0), = (0, v7), U = p = 
W = ^/J /r, and v = M/ir'^y/J), where we use the tilde to avoid mixing of the variables inside the star. From the 
junction conditions at the surface of the dust sphere, we obtain the relationships between the (77, x)-coordinates and 
the {t, r)-coordinates, given by 



Ts = i?(?7)sinxo, 



t 

2M 



= In 



[(r,o/2M)-l]i/2 + tan(r7/2) 



[(r«o/2M) - 1]V2 - tan(j?/2) 



\2M 



1/2 



f— ) 
V4M/ 



{rj + sin 77) 



(2.35) 

(2.36) 



where Tso = rs{t = 0) = 2M/ sin^ xo is the initial stellar radius in Schwarzschild coordinates. 



C. Magnetic Fields 

As mentioned the above, we consider weakly magnetized dust spheres in which the magnetic effects are treated as 
small perturbations on the OS solution, where we can consider that the electromagnetic fields are axisymmctric due 
to the nature of a spherical symmetric background. Thus, perturbations of the electromagnetic fields, dF/^,^, can be 
described in terms of the spherical harmonics F'"* by the following relations 

6F01 - -SFio = 62^'" , (2.37) 

6Foa = -5F,o = eiS'Z + es^^'™ , (2.38) 

SF,a = -6Fai = hS'Z + hY'Ta , (2-39) 

5F23 = -SF32 = 62£23i^'" , (2.40) 
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where S' ™ = £ b- These variables of are governed by the Maxwell equations, i.e. 



(2.41) 
(2.42) 



where SJ^ is the perturbations of the current four-vector. The equations (|2.4ip and ()2.42p are correct up to order of 



In the interior of the dust sphere, we consider infinitely conductive fluids, i.e., the ideal magnetohydrodynamic 
approximation has been adopted, according to which SF^^u'^ = 0. With this approximation and Maxwell equation 
(|2.4ip . we obtain the basic equations for electromagnetic fields inside the star such as 



ei = 62 = 63 = . 



l{l 



l)bi+d^b2 = 0. 



(2.43) 
(2.44) 
(2.45) 



The equation (|2.44p tells us that the magnetic field does not change for a comoving observer. Additionally, with the 
Faraday tensor F^"^, since the magnetic field, S^, can be described as = Cp^uapu'^ F"^ one can see that the 
quantities 6i and 62 are related to the poloidal magnetic component while 63 is associated with the toroidal magnetic 
component. Furthermore, as similar to Paper I, in this paper we consider the case that the magnetic fields are confined 
in the stellar interior, i.e., 61 = 62 = 63 = at stellar surface. (In fact the necessary condition that the magnetic 
fields are confined inside the star is only 62 = at stellar surface, which is derived from the junction conditions 
for the magnetic fields as Paper I, but for simplicity in this paper we adopt the above conditions.) Of course the 
electromagnetic fields outside the star are also important to see the emission of electro-magnetic waves, but they can 
be seen in the near future elsewhere. Additionally in this paper we focus only on dipole electromagnetic fields, i.e., 
electromagnetic fields associated with I = \. In the next section we will see how the dipole electromagnetic fields can 
drive the quadrupole gravitational radiation, which is more important in the observation. 



III. PERTURBATION EQUATIONS FOR POLAR PARITY 
A. Interior region of the star 

The perturbation of the energy-momentum tensor St^^^, is described as 



(3.1) 



where 5t^^^ and 5t^^J^'^ are corresponding to the energy-momentum tensor for the dust and the electromagnetic field, 
respectively. Since the polar perturbation of fiuid 4- velocity, 5u^, and the perturbation of density, 5p, are defined as 

r / 1 

5u„, 



the expansion coefficients of St^/^^ in Eq. (|2.5p are 



^^AB — P 



7 {uAnu + tiaub) + ^ {^bcua + hAcUB) 



LUpUAUB, 



3(M) 



(3.2) 
(3.3) 

(3.4) 

(3.5) 
(3.6) 



On the other hand, as mentioned before we consider the effect of the dipole magnetic fields on the quadrupole 
gravitational radiation. In this case the non-zero expansion coefficients for St^i^'^ associated with I = 2 gravitational 
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waves are given in the Appendix A of Paper I as 



A/(EM) 



Ai2(EM) 
^^3(EM) 



bi' 



1^ 
b2 



R? sin X 

2 



sin X 



ii^sin X. 
sin^ X 



siv? X ' 
n (53' - 



i?2 



K&2 



(3.7) 
(3.8) 
(3.9) 

(3.10) 
(3.11) 



i?4 sin-^x' 

where k = (SttVStt)^^. Thus we can derive the gauge- invariant quantities for the total matter perturbations such as 



Tip, = ~phpi + i^pR 



T — 



b2' 



bi 



R^sin X sin X 



bi 



R^ sin X sin x 



T 
T 

T 
^x 



-jpR^ 
—apR, 

R^ siv? X ' 

At 62' 



i?4 sin^x' 

where a, 7, and w are gauge-invariant set of fluid perturbation defined as 



A 

a — a ~ p UA, 

A 

7 = 7 — n 



1 



P^u^lB + [Pb\a -Pa\b) 



Cb-p (Inp)i^. 



(3.12) 

(3.13) 

(3.14) 
(3.15) 
(3.16) 

(3.17) 
(3.18) 

(3.19) 
(3.20) 
(3.21) 



Then, with Eqs. (|2.19|) - (|2.21[) . the evolutionary equations for the metric perturbations on the interior region are 
described as 



R 
R "^"^ 



smx 



cosx 
sinx 



- 2fc - 



sin^ X 



2 sin X 



(3.22) 
(3.23) 
(3.24) 



where we set to be Cg = since the matter is considered as dust. Notice that the source terms Sq, Sk, and are 
produced only by the perturbed magnetic fields and the concrete forms are given in Appendix [X] In other words, as 
mentioned in ^11 A 21 we can calculate the metric evolutions apart from the matter perturbations. At last, with Eqs. 
(|2.23p - (|2.25p , the evolutions of matter perturbations are determined by using the variables for metric and magnetic 
perturbations, which are described in Appendix IB] From the above system of equations for metric perturbations, we 
can see that the variable C is independent of the others k and ip, while the equations for variables k and ip do not have 
such terms as d-^k and d^ip. Thus for the interior region it is enough to calculate the evolution for only After that. 
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in order to adopt the junction conditions (see ^III C[) . we have to calculate k and tp in the vicinity of stellar surface 
with the determined C- 

Considering the behavior of the metric perturbations near the stellar center, we introduce a new variable, C, which is 
regular at the stellar center and defined as C = (-Rsinx)'^^ ^ [IHl- With new variable the above perturbation equation 
for (, Eq. (I3.22p . can be rewritten as 



^d'^C + dlQ - 2{i + 3)^a,c + 2{i + i)^a,c - + 2) 

' * R smx 



dr,R 



+ 1 



c = 



Sr 



(Rsinx) 



1+2 ■ 



(3.25) 



Furthermore, in the actual numerical calculations we adopt the double null coordinates, (u, v), defined as u = rj — x 
and V = rj + X- In these coordinates, the perturbation equation is rewritten as 



dudv 



R 



smx 
cos X 



(/ + 3)^-(Z + l) 

R smx 



du 

ac + 2) 

dv 4 



d^R 



+ 1 



Sc 



4 (i?sinx) 



1+2 ■ 



(3.26) 



B. Exterior region of the star 



In the exterior region the master equation for perturbations can be reduced to the well known Zerilli equation for 
the Zerilli function, Z(t,r), which has the form 



Vz{r) = / 



-dtZ + dl^ Z - Vz{r)Z = 0, 
1(1 + 1) 6M r^X{X + 2) + 3M{r- M) 
(rA + 3M)2 



(3.27) 
(3.28) 



where X = {I + 2)(l — l)/2 and the tortoise coordinate, is defined as r» = r + 2M \n{r /2M — 1), which leads to 
dr, = fdr- This variable Z{t,r) is constructed with the variables ( and k as 



where 



Z{t, r) ^ A{r)C{t, r) + B{r)k{t, r) + C{r)dr, k{t, r). 



2r rAi -2r^ 

A{r)^-, B{r)^-^, C{r)^-— 
A3 A3 /A3 



(3.29) 



(3.30) 



where Ai = -1 + {P + I + l)/f and A3 = -3 + (P + / + 1)// '2^. On the other hand, the variables of metric 
perturbations, fc, C and 'ip, are produced with Z as 



k^^ 



r l{l + l) 



6MZ 



dr,Z 



(3.31) 



~_ 4r + 3M 
C — — ^r^rrzZ 



QM 4r2 ~ 7rM ~ 9AP 



2M 



5r - 3M 



r(2r + 3M)" ' ;(? + l) r^2r + 3M)^ " l{l + 1) {r - 2M){2r + SM)"''''" ' l{l + 1)/^'''^' ^^'^^^ 



-dr.Z 



2r 



1(1 + I) f 



{l-l){l + 2) M\ ■ 

7 > OtZ + rdtdr,Z 

A3 r J 



(3.33) 



Finally, the Zerilli equation (|3.27p can be rewritten in terms of the double null coordinates, u = i — r* and -C = t + r*, 

as 



(3.34) 



C. Junction Conditions at the stellar surface 



In order to connect the metric perturbations of interior region with those of exterior regions, we have to impose 
the junction conditions at the stellar surface, which are derived from the continuity of the induced metric and the 
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extrinsic curvature [261]. As mentioned the above, in this paper we focus on the case that the magnetic field would be 
confined inside of star, i.e., q — a,t stellar surface (see Ea. (|A7p ). In this case the junction conditions to impose are 
the continuity of N, fc, ^7 "0, k' + SirpN, and (' + 2/^-0, where the variable N is determined by using 7 and such as 
N — = — (7 + -0/2). These conditions are equivalent to 

(^(^ k^k, Tp^^p, (3.35) 

OR / \ Fl / \ 

d^k + 87:pRN = sinxo [dtk) + j cosxo [dr,~kj , (3.36) 

dxC^^ sin xo {dtC) + J cos xo {dr, c) ■ (3.37) 

With these junction conditions, the exterior variable Z can be described by using the interior variables ^, fc, and tp 
I26il. such as 



2r4 1 

^ — j,k-\ 

(; + 2)(Z- l)r + 6Mi?2 




(C + fc) + _ (a^fc + 8.pi?7V) + 

smxo -fl smxo -fi 

(3.38) 



IV. NUMERICAL PROCEDURE 

The numerical procedure adopted in this paper is basically similar to that in Paper I. In order to calculate the 
nonspherical perturbations, we divide the background spacetime into three regions named I, II, and III (see Fig. 
H]). Region I denotes the stellar interior while regions II and III are corresponding to the exterior. Region II is 
corresponding to the intermediate exterior region, which is introduced to help the matching procedure at the stellar 
surface in numerical computation. Region III is separated from region 11 via the null hypersurface defined by w = wq, 
which is the ingoing null ray emitted from the point where the stellar surface reaches the event horizon, i.e., the point 
1-L in Fig. m 

In order to solve the wave equations numerically, we adopt the finite difference scheme proposed by Hamade and 
Stewart [3^], in which we use the double null coordinates (u, v) in region 1 and (m, v) in regions II and 111. In region I, 
to avoid numerical instabilities we integrate the wave equations by using a first order finite difference scheme, while 
in regions II and III the numerical integration is a second order finite difference scheme. In region 1 we adopt the 
equally spaced grids for (u,w), i.e.. Aw and Aw are constant, and we set to be Am = Aw. With this assumption, the 
interval for 77 and X) A77 and Ax, are also constant and A77 ~ Ax. The grid points in region 11 are determined so that 
at the stellar surface they agree with the grid points produced with the coordinates in region I. That is, Am and Au 
are not constant and Am ^ Av in region II. It is noted that with initial data sets on 77 = in region I and on m = mq 
in region II, the evolutions in regions I and II can be calculated, independently of the information in region III. After 
the calculation of the evolutions in regions I and 11, with the data set on the null hypersurface, w = wo, and initial 
data set on t2 = mq the evolution in region 111 can be calculated, where we adopt the equally spaced grids for (m,^). 
As mentioned before, hereafter we focus only on the quadrupole gravitational waves (Z = 2), which are coupled with 
the dipole magnetic fields {I = 1). 

Finally it should note that in region 111 we calculate the time evolution for only variable Z , which is subject to the 
Zerilli equation (|3.34p . while in region 11, to make it easy to deal with the junction conditions on the stellar surface, 
we also calculate the time evolutions for the variable C as well as Z . The perturbation equation for C, is derived from 
the equation (|2.19p . which is described with null coordinates as 



9\ 
dudv 




(4.1) 



In the rest of this section, we describe the initial data, the boundary conditions at the stellar center and at the spatial 
infinity, and the special treatment of the junction condition as the stellar surface approaches the event horizon. 



A. Initial Data 

To start numerical simulations, we need to provide a data set on the initial hypersurface for the quantities C,, duC, 
and dyC as well as the magnetic perturbations of &2 and 63 for the interior region, while Z, duZ , and dyZ for the 
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f] = n 

stellar surface 



FIG. 1: A schematic description of the Oppenheimer-Snyder spacetime for the collapsing model in characteristic coordinates. 
Region I denotes the stellar interior while regions II and III correspond to the exterior. The stellar surface, where r = or 
X ~ Xo, is the boundary between regions I and II, and the shaded region corresponds to the stationary region outside star. 

exterior region. Outside the star, we assume that the initial perturbations are "momentarily static" , which is similar 
initial condition in [13, [131 . With this assumption the initial distribution of Z{r;t = 0) is determined by using the 
following equation; 



dlZ + Vzir)Z = 0, 



(4.2) 



with the boundary condition at infinity as 



Z{r;t = 0) 



92 



'2M 



(4.3) 



where (72 is a constant denoted the quadrupole moment of the star. Similar to |20l . \2l\ . we assume that (72 = 2M. 
Since this solution is a static, the initial perturbation outside the star, Z{r), does not evolve until a light signal from 
the stellar interior arrives there, i.e., on the gray region in Fig. [T]the solution Z{r) will not be changed. Thus we 
can use the initial data Z{r) as the data set on the null hypersurface u = uq. Furthermore, with the assumption that 
dtZ = at i = 0, the data for duZ and d^Z are given as duZ = —{dr,Z)/2 and dyZ = {dr,Z)/2, respectively. 

With respect to the initial condition inside the star, we can choose appropriate functions of magnetic distributions, 
62 (x) and &3(x)j where the conditions to determine the electromagnetic perturbations are Eqs. p.43p - (|2.45|) . As 
mentioned before, since in this paper we focus only on the case that the magnetic fields are confined inside the star, 
we should put the boundary conditions at x = Xoi such as &i = &2 = ^3 = 0. Then, similar to the exterior region, if 
the momentarily static condition for C would be assumed, with the given initial distributions for magnetic fields, the 
initial data for C can be determined by integrating the equation of 



die 



6 cosx 
sinx 



(i?sinx) 



(4.4) 



It should notice that for the case of the non-magnetized sphere the value of d^C,/C, at the stellar surface is independent 
from the central value of because the equation (|4.4[) does not have the source term. Thus in this case we produce 
the initial data of C, so that at the stellar surface the metric perturbation is not smooth but just continuous. However, 
the effect of this non-smoothness on the emitted gravitational waves looks like very small (see Figs. [5] and [3). 
Actually Cunningham, Price & Moncrief also adopted the non-smoothness initial condition in their calculations |20j. 
Additionally the initial data for fc, dy^k, and at x = Xo are derived from the junction conditions. Then, as mentioned 
in the previous section, the variables for matter perturbations, 7(x)i ol{x), and lij{x)i can be determined by using the 
initial distributions for metric perturbations. Finally we also add an assumption that TV = at 77 = 0. 
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B. Boundary Conditions 



For the numerical integration we have to impose the boundary conditions. One is the regularity condition at the 
stellar center (x = 0) and the other is the no incoming-waves condition at the infinity. The regularity condition at 
the stellar center demands that d^^C, = 0, which is reduced to d^C = dyC,. With respect to the no incoming radiation 
condition at the infinity, we adopt the condition as dZ/du — (see, e.g., [33|). 



C. Special Treatment of the Junction Conditions near the Event Horizon 

When the stellar surface reaches the event horizon, the junction conditions discussed earlier in §111 CI can not be 
used any more because the terms related to /^^ diverge. Instead of these junction conditions, following [2^, we adopt 
an extrapolation for the value of Z on the junction null surface, v = vq, in the vicinity of the point Ti in Fig. [T]as 

^ = ^""'"+ >„X. (r^r"-), (4.5) 

. Z--- + {r-- - , (4.6) 

where Z" and r" are the values of Z and r on v = vo at n-th time steps, while A'max denotes the total number of 
time steps in region II, and r™ — 2M. 



V. CODE TESTS 



In order to verify our numerical code, we have calculated quadrupole gravitational radiations emitted during the 
collapse of a nonmagnetized homogeneous dust sphere, i.e., the gravitational waves emitted from the Oppenheimer- 
Snyder solution. The number of spatial grid points inside the star, N^, which corresponds to region I, is chosen 
to be N-^ = 1000, because we can not see a dramatic improvement with lager number of grid points. Actually, as 
shown in Fig. ^ the waveforms of gravitational waves emitted from the collapsing dust ball with Ny, — 1000 are very 
similar to those with = 2000, and the total energies of emitted gravitational waves defined later also agree with 
each other within 0.626 % for r^o — 8M and 1.37 % for r^o = 20M. In region III, the step size for integration is 
given as Au = (wmax — uo)/Nu, where Umax is determined with the expected maximum time for observer, tmax, and 
the position of observer described in tortoise coordinate, r*ob, as Umax = ^max — ''*ob- In this paper we adopt that 
^max = 2000M and rob = rso + 40M, respectively, where the position of observer is the same choice as the previous 
study by Cunningham, Price & Moncrief 20] . Since the numerical code in this region is essentially the same as that 
in Paper I, the number of grid points for outgoing null coordinate, u, in region III is assumed to be Nu = 10000 in this 
paper (see Table I in Paper I for the convergence test). Then we have only one parameter to determine the emitted 
gravitational waves, i.e., the initial radius r^o. 

As noticed in [20[, the emitted gravitational waves are characterized by the quasi-normal ringing oscillation and 
subsequent power-law tail. In Fig. [31 we show the waveform of the gravitational waves for / — 2 emitted during the 
collapse of the homogeneous dust ball, where the left and right panels are focused on the quasi-normal ringing and 
on the power-law tail, respectively. The fundamental frequency of the quasi-normal ringing have been calculated by 
Chandrasekhar & Detweiler [s^, such as 2Muj = 0.74734 -I- 0.17792i. On the other hand, our numerical results show 
that the oscillation frequency is 2Muj = 0.737, which agrees well with the previous value with only 1.3% error, while 
the damping rate also consorts with the theoretical value (see the left panel of Fig. [3]). As for the late-time tail, 
in the right panel of Fig. [31 we find that the amplitude of gravitational wave decays as {t — io)~^: where Iq is the 
time when the observer receives the first signal emitted from the stellar surface, i.e., Iq = r»ob ~ '''*sO- This result 
is in good agreement with the analytical estimate by Price [s^, that is {t — to)~'-'^''~^^^ ■ Through these estimations 
for the frequency of quasi-ringing and the late-time tail, we believe that our numerical code is possible to derive the 
gravitational waves with high accuracy. 

At the end in this section, we compare the total energy emitted during the collapse with the previous results by 
Cunningham, Price, & Moncrief [13] (CPM1979). It is worth to notice that the variables inside the star adopted in 
CPM1979 axe different from those in the equation system for the gauge-invariant formalism proposed by Gerlach & 
Sengupta f22] and Gundlach & Martm-Garcia [26[. The total emitted energy, -EgWi is estimated by integrating the 
luminosity of gravitational waves, icw, with respect to time, where the luminosity is defined as 

Lgw = :^{Z,)\ (5.1) 
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FIG. 2: Waveforms of gravitational waves for I = 2 emitted from the homogeneous dust collapse with the initial radius r^o ~ 8M 
(left panel) and r^o ~ 20M (right panel). The solid and broken lines are corresponding to the results for A^^^ = 1000 and 2000, 
respectively, where A''^ is the number of spatial grid points inside the dust ball. 
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FIG. 3: Waveforms of the quadrupole gravitational radiation emitted during the collapse of the non-magnetized homogeneous 
dust, as function of time. The initial radius of the dust sphere is set to Vso = 8M while the fiducial observer is set at 
Tob = fso + 40M. In the right panel the late-time is compared with its theoretical value {t — to)~'^'^'^' , where to is the time 
when the observer receives the first signal emitted from the stellar surface. 



for the I — 2 gravitational waves (e.g., [1^). In Fig. [4] we show the total emitted energy of gravitational waves as a 
function of the initial stellar radius Tso, where for comparison we also plot the result of CPM1979. It notes that in 
this figure we adopt the normalization for the quadrupole moment, q2, so that q2 = 2M as mentioned before. This 
figure shows that there are small difference between our results and those obtained in CPM1979. The main reason for 
this difference could be the difference how to choose the variables inside the star. Additionally as we noticed in Paper 
I, the difference of the accuracy in numerical code might be also added in the reason. Anyway we can observe that 
the total emitted energy systematically decreases as the initial stellar radius increases and that the emitted energy is 
very similar to that of CPM1979. Furthermore, with our variables inside the star and with our initial data, we derive 
the empirical formula for the total emitted energy as 

^=3.37xl0-^x f^r', (5.2) 
2M \2M) ' ^ ' 



which is also plotted in Fig. [4] 
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FIG. 4: Total energies emitted in gravitational waves from the homogeneous dust collapse without magnetic field as a function 
of the initial stellar radius, where rso = 6M, 8M, 12M, 16M, 20M, and 24M. The filled circles correspond to the results by 
Cunningham, Price, and Moncrief [20[ |. while the squares are corresponding to our numerical results. The broken line denotes 
the empirical formula derived from our results. 



VI. GRAVITATIONAL RADIATIONS FROM THE OS SOLUTION 



In order to calculate the gravitational waves emitted from the collapsing phase of a magnetized dust sphere, we 
have to provide the initial distribution of the magnetic field. In other words, one needs to set up the functional forms 
of 62 and 63 on the hypersurface rj = 0. The initial distributions can be determined as the following two conditions 
are satisfied; (a) the regularity condition at the stellar center and (b) the junction condition at the stellar surface. 
Since we made assumptions in this paper that the magnetic field is confined inside the star and then the value of q 
becomes zero at stellar surface, the conditions at the stellar surface can be described as 61 (xo) — b2{xo) — bsixo) ~ 0. 
Now we introduce two new variables, 62 and 63, such as 



hix) = S2sin^x^2(x), 
hix) = S3sin3x^3(x), 



(6.1) 
(6.2) 



where B2 B3 are arbitrary constants related to the strength of the magnetic field. With analytic functions 62 and 63, 
the regularity condition at the stellar center for the magnetic field is automatically satisfied. Since the geometry of 
the magnetic field when the collapse sets in is practically unknown, in this paper we adopt the following two types of 
the initial distributions for the magnetic field; 



il):b2{x) = h{x) = l-2 



(II) : 62 (x) = bsix) = 16{ — 

\Xo 



1-2 ( ^) +(^^ 



Xo 



Xo 



(6.3) 
(6.4) 



where the maximum value of &2 and 63 are chosen to be one in the range of < x ^ Xo- For the first profile (I) the 
magnetic field is stronger in the center of the sphere, while for the second profile (II) the field becomes stronger in the 
outer region. Additionally it notes that for both profiles the value of bi{x), defined as 6i(x) — ~f^x^2(x)/2, becomes 
zero at the stellar surface, i.e., as mentioned before the value of q is zero at x = Xo- With these magnetic profiles, we 
found that the allowed values for B2 and B3 have to be in a part of the range of B2 < B3, in order to produce the 
initial data set so that the inner metric perturbation should be smoothly connected to the stationary solution in the 
outer region at the stellar surface. So, in what follows, we consider the two cases for magnetic field; one is that only 
toroidal magnetic component exists, i.e., B2 = 0, and second is that the poloidal magnetic component also exists as 
well as toroidal one, where those are satisfied the condition that B2 < B3. 
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A. Toroidal Magnetic Field 



First we consider the case that only toroidal magnetic component exists, i.e., B2 = 0. In this case, the source term 
in the equation (j4.4p to determine the initial distribution, ^, is proportional to B^^. The value of S3 is determined 
so that the initial inner metric perturbation should be smoothly connected at the stellar surface to the stationary 
solution for exterior region. Then we can get the distributions for the initial inner metric perturbation, Cix)i ^^nd 
for the initial density perturbation, which are shown in Fig. [5] with the two different magnetic profiles (I) 

and (II). In this figure the initial stellar radius is set to be Tsq = 8M, but the functional forms of i^(x) and w(x) 
with the different initial stellar radii are very similar to that with r^o — 8M. From this figure, we can see that the 
initial distributions of ( and io depend strongly on the magnetic profiles even if the initial metric perturbations for 
exterior region are adopted the same as the stationary solution with (72 = 2AI. On the other hand. Fig. [6] shows 
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FIG. 5: The initial distribution of inner metric perturbation, C(x)) on the left panel and the initial density perturbation, a;(x), 
on the right panel, where the initial radius is adopted that Vso ~ 8M. 

the waveforms of the emitted gravitational waves with these initial perturbations, where the left and right panels are 
corresponding to the results with the initial radius r^o = 8M and 20M, respectively. For comparison, the waveforms 
for the non-magnetized dust collapse are also plotted. The first observation of Fig. [6] is that the waveforms of the 
emitted gravitational waves are almost independent from the magnetic profiles in spite of the difference of initial 
perturbations inside the star. That is, if the initial dust sphere consists only of the toroidal magnetic component, 
it might be difficult to distinguish the interior magnetic profile by using the direct detection of the waveform of the 
emitted gravitational waves. Additionally we can observe the difference between the waveforms of gravitational waves 
with the toroidal magnetic field and without magnetic field. With smaller initial radius, the shape of waveform is 
similar to that for the nonmagnetized case, still we can see the effect of the existence of magnetic field, i.e., the quasi 
normal ringing can be seen earlier and the amplitude is also enhanced a little due to the magnetic effect. While, 
with large initial radius, it is possible to watch the obvious influence of magnetic field on the waveform of emitted 
gravitational waves, where the amplitude of waveform grows large and the maximum value of gravitational wave 
becomes a negative. In other words, with large initial radius, the waveform before the quasi normal ringing would 
be observed can be changed remarkably. The reason for this could be thought that with large initial radius it takes 
longer time until the stellar surface reaches to the event horizon and then the inner magnetic field can affect on the 
metric perturbations with longer time. In particular, the dependence of magnetic effect on the initial radius can be 
seen clearly in the total energy of emitted gravitational waves. Fig. [7] shows the total energy as a function of the 
initial stellar radius with the circles for magnetic profile (I) and with the triangles for magnetic profile (II) , where for 
comparison the total energies for the nonmagnetized dust case are also shown with the squares. It is found from this 
figure that with large initial radius, due to the magnetic effect the total energy of emitted gravitational wave becomes 
much larger than that for the nonmagnetized dust collapse. Then the dependence of the total energy for the dust 
collapse with toroidal magnetic field on the initial stellar radius, is quite different from the empirical formula (|5.2p for 
the nonmagnetized dust collapse. 
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FIG. 6: Waveforms of gravitational waves for I = 2 emitted from the homogeneous magnetized dust collapse with the initial 
radius rso = 8M on the left panel and with Vso = 20M on the right panel. The solid and dotted lines are corresponding to the 
results with the magnetic profile (I) and (II), respectively, while for comparison the result for the non-magnetized case is also 
denoted with the dot-dash line. 
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FIG. 7: Total energies emitted in gravitational waves from the homogeneous magnetized dust collapse as a function of the initial 
stellar radius, where Vso = 6M, 8M, 12M, 16M, 20M, and 24M. The filled circles and triangles correspond to the calculated 
results with magnetic profile (I) and (II), respectively, while the squares are corresponding to the total energy emitted from 
the non-magnetized dust collapse. 



B. Poloidal and Toroidal Magnetic Fields 



Next we consider the magnetic field, which consists of the poloidal and toroidal components. In this case we can 
introduce the new parameter, f}, defined as /? = B2/B3, and if we choose the value of (3 the initial inner metric 
perturbations are determined as those should be smoothly connected to the outer stationary solution. It notes that 
the case for /3 = corresponds to the dust model, which only toroidal magnetic component exists shown in the 
previous subsection. Table U shows the allowed maximum values of /3 with different combinations of magnetic profiles 
for the poloidal and toroidal components and with different initial stellar radii. From this table, it can be seen that 
with large initial radius it becomes more difficult to produce a magnetized dust model with large value of and that 
the maximum values of f3 depends strongly on the inner magnetic profiles. 

Similar to the collapse of magnetized dust with only toroidal component, the magnetic effects can be seen in the 
waveforms of gravitational waves a little. Fig. [5] shows waveforms for the collapse of magnetized dust with = 8M 
and with several values of /?, where the magnetic profile (I) is adopted for the poloidal and toroidal components. From 
this figure it is found that the emitted gravitational waves are basically characterized by the quasi-normal ringing as 
well as the case of nonmagnetized dust collapse. While, we can also see the specific magnetic effects in the waveforms, 
where as the value of (3 becomes larger, the amplitude of gravitational waves is enhanced and its maximum value 
changes from a positive to a negative. These are similar features to the case of dust collapse with only toroidal 
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TABLE I: Allowed maximum values of /3 with the diflterent combinations of magnetic profiles for the poloidal and toroidal 
components and with different initial stellar radii. As mentioned before, 62 and 63 are associated with the poloidal and toroidal 
magnetic components, respectively. 



profile 


TsO 


62 63 


6M 


8M 


12M 


16M 


20M 


24M 


(I) (I) 


0.073 


0.054 


0.035 


0.026 


0.021 


0.017 


(I) (11) 


0.26 


0.20 


0.13 


0.10 


0.080 


0.067 


(11) (I) 


0.020 


0.014 


0.0094 


0.0070 


0.0055 


0.0046 


(11) (II) 


0.072 


0.053 


0.035 


0.026 


0.021 


0.017 



component with large initial radius. In other words, with large value of magnetic ratio, even with small initial radius, 
we can see the magnetic effect in the waveform before the quasi-normal ringing would be observed. This tendency 
holds for the magnetized dust collapse with different magnetic profiles. As a result, with large value of magnetic 
ratio, the total energy of gravitational waves grows. The total energies for rso = 8M and 20M are plotted in Fig. [9] 
as a function of magnetic ratio, where the different lines are corresponding to the different combinations of magnetic 
profiles for the poloidal and toroidal components. In the explanatory notes of this figure, for example, "(I)-(II)" 
denotes that the magnetic profiles (I) and (II) are adopted for the magnetic variables, 62 and 63, respectively. This 
figure tells us that the total energies emitted in gravitational waves from the magnetized dust collapse depend strongly 
on the magnetic ratio and the inner magnetic profiles. This sensitivity, as well as the change of waveforms due to the 
existence of magnetic fields, could be important to extract some information of inner magnetic profiles of progenitor 
from the direct observation of gravitational waves during the black hole formation after the stellar collapse. 
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FIG. 8: Waveforms of gravitational waves emitted from the homogeneous magnetized dust collapse with the initial radius 
Vso ~ 8,M on the left panel, where three different lines are corresponding to the dust models with different magnetic ratio, /3, 
and the magnetic profile (I) is adopted for both poloidal and toroidal components. It notes that the line for /3 = is the result 
of collapse without poloidal magnetic component. The right panel is the magnification of the region encompassed by the square 
in the left panel. 



VII. CONCLUSION 



In this article, with the gauge-invariant perturbation theory we have studied the dependence of the stellar magnetic 
fields on the polar gravitational waves during the collapse of a homogeneous dust sphere. It should be emphasized 
that this is first calculation of emitted polar gravitational waves on the dynamical background spacetime with the 
covariant gauge-invariant formalism on the spherically symmetric spacetime and the coordinate-independent matching 
conditions at stellar surface, which is devised by Gundlach and Martm-Garci'a [2^. So far, such calculations could 
not be done due to the difficulty to treat the boundary conditions at stellar surface. While, in order to solve this 
difficulty, we evolve not only Zerille function, Z, but also metric perturbation, in the intermediate exterior region 
(region 11) and the calculation of emitted gravitational waves is successful. 
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FIG. 9: The total energies emitted in gravitational waves from the homogeneous dust collapse with magnetic field as a function 
of the magnetic ratio, /?, where the left and right panels correspond to the results for the initial radius Vso = 8M and for 
Tso = 20M. In the figure, the different lines are corresponding to the different combinations of magnetic profiles for the poloidal 
and toroidal components. In the explanatory notes, for example, "(I)-(II)" denotes that the magnetic profiles (I) and (II) are 
adopted for the magnetic variables, &2 and 63, respectively. 



With this nemerical code, we consider the magnetic effects on the polar gravitational waves from the Oppenheimer- 
Snyder solution describing collapsing dust, where the magnetic fields are introduced as a second order perturbation 
term. Even if the initial magnetic perturbations are small, as the collapse proceeds they could get amplified and 
become significant because of the conservation of magnetic flux. In particular, similar to Paper I, we have assumed 
that the magnetic field is axisymmetric, where the dipole magnetic field perturbations are the ones that couple to the 
quadrupole polar perturbations of the gravitational field. Additionally, we assumed momentarily static initial data 
and we have not taken into account the influence of the exterior magnetic field in the propagating gravitational waves. 

Through the investigation, it is found that there exists an evidence of the strong influence of the magnetic field in 
the gravitational wave luminosity during the collapse. Depending on the initial profile of the magnetic field and its 
ratio between the poloidal and toroidal components, the energy outcome can be easily up to a few order higher than 
what we get from the nonmagnetized collapse. In addition, it is possible to observe an important change before the 
quasi-normal ringing is detected, which is induced by the presence of the magnetic field. These magnetic effects can 
be seen the collapsing model with large inital radius and with large magnetic ratio between the poloidal and toroidal 
components, since for a large initial radius the time needed for the black hole formation is longer and then the 
magnetic field acts for longer time on the collapsing fluid. It notes that the magnetic effects on the polar gravitational 
waves are different from those on the axial ones, i.e., the axial gravitational waves are independent from the magnetic 
ratio while they depend only on the magnetic strenght such as the value of x S3. Such magnetic effects could be 
helpful to extract some information of inner magnetic profiles of progenitor form the detection of gravitational waves 
radiated from the black hole formation after the stellar collapse. 

At the end, we believe that although this study might be considered as a "toy problem" it has most of the ingredients 
needed in emphasizing the importance of the magnetic fields in the study of the gravitational wave output during 
the collapse. The final answer to the questions raised here would be provided by the 3D numerical MHD codes 
(see in [isl |37| for the recent developments), but this work provides hints and raises issues that need to be studied. 
Furthermore, as future studies, we condiser to study magnetic effects on the gravitational waves emitted from the 
more complicated background collapsing models such as inhomogeneous dust collapse and the stellar collapse with 
perfect fluid, while it should be also important to take into account the background magnetic field such as [38| . 
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APPENDIX A: CONCRETE EXPRESSION FOR THE SOURCE TERMS 



In this appendix, we show the concrete expression for the source terms in the perturbation equations, which are 
not written in the main text. Those for Eqs. - (P?^ and (P?^ - (1^?^ are 

5c =167r [T^ + 2{n^TAy - n^n^TAB + 2(2i/ - W)n'^TA\ + 4:{fi - U)k + 3K + - bv)C' 
- 2{fi - U)q + (8i^ - &W)q' + 2q" + 2 [2(/i - U)W - 2fiiy - n' + i> + STrt^-^w^ns] V 
2 



2W' - 2iy^ 



4W' + iiyW - 2W^ - 2^ 



/(/ + !) + 4 
2r2 



167rQ + 87r(i'^^ +t'4^PAi3) 



(Al) 



5fc =87r [(-c^u^ii^ + n^n^)TAB + mn^TA] 

+ (c^ + 1)[/C + [4C/ + cl{^i + 2[/)] fc + (c^ - l)WC' + 2clW)k' - 2Uq + 2Wq' 



%TT{t\ ~ ^'"pab) + 2U{2^i + U) + 



8n{t\ + f'^PAB) + 2[W'-^ 



iC + k) 



a-i)a + 2) 

2^2 



{8TT{t\ - t^^PAB) + 2U{2fi + U)} - 8n{t\ + ^^^pab) + &W 



+ 2 [{cl + l)U{v + W) + (c2 - l)^iW - 8tt{cI + l)i^^w^7is] 
^2v{C + fc) + 2^1/^ + C' - 2q' + 2(M^ - - IQim'^TA, 



(A2) 
(A3) 



and 



=?/C + (m + 2C/)fc + T4^C' - 2iyfc' 

?7r(t^^ - t^^PAij) + 2J7(2^ + U) + 



/(/ + !) + 2 
2r2 



+ C/(2/i + C/) - VF(2t/ + W^) + STTt 



[87r(t^^ - t^^pAB) + 2C/(2/i + [/)] q + 2 [i^C/ + ^iW + UW - 8TTt^^UAnB] ^ 



Ca =2^(C + k) + 2viP + C + 2k- 2q{fi + U), 
where q is given by Eq. (|2.22p . such as 



(A4) 



(A5) 
(A6) 



J, 2 7 2 

oi - 03 



(A7) 



Additionally, the source terms in the perturbation equations for the metric perturbations p.22p - (|3.24p are 



Q 6cosx„ , 4cos2x-^(? + l) + 4 
Sr = 2d^q ■. d^q + ^ q 
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hi 
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sinx 



c o;^ 2 cosx 
= 20yq : q - 
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26162 



57r_R2 sin x 



(AS) 
(A9) 

(AlO) 
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APPENDIX B: MATTER PERTURBATIONS FOR INTERIOR REGION 

The matter perturbations are given by using the variables for metric and magnetic perturbations, which are deter- 
mined after calculation for evolution of metric perturbations. With Eqs. (|2.23p - (|2.25p . those are 



1 



smx ^ R R ^ R |2sin2x 



SirpR^ 



* R R R smx 



• 2 

sm X 



2 sin X 



R 



2cosx 
sinx 



d^R 



1 - AirpR^ } ^ 



d^k 



(Bl) 



AdjjR cos X 



R sin X i?2 g jj^^ ^ 



sin^X 



1 



WnpR 



+ ^ (C + fc) + + 2a,fc ^^"^ ' 



R 



R 



(B2) 
(B3) 
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